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1. Introduction 

In this paper we wish to study the conditions for isotropisation of the Generalized Scalar- 
Tensor theory plus a massive scalar field in the Bianchi type / model. They are many 
reasons to be interesting by this model 

First of all General Relativity is a good description for the weak gravitational fields 
(solar system tests) as for strong ones (binary pulsar) although deviating are expected 
for early times or in extreme cases such as black hole. Then, it is interesting to consider 
a Lagrangian whose "geometric" part looks like General Relativity. Moreover particle 
physics progress and the idea of infiation in the eighties show that scalar fields could 
be essential components of a gravitational theory. In this work we will consider a 
massive one. It could be justified by observations of type I A supernovae [H [2| which 
seem to demonstrate that the dynamical behaviour of our Universe is accelerated. 
Most of time this is interpreted as the presence of a cosmological constant in the field 
equations although other explanations can be advanced such as this of a non-perfect 
fluid with quintessential matter [H]. The Boomerang experiment jl| indicates that it 
could represent the dominant energy of our Universe. However, the present value of this 
constant is in disagreement with this predicted by particle physics at early times. One 
way to solve this problem is to consider a varying potential U, i.e. a massive scalar field. 
We will also describe the coupling of the scalar field with the metric functions by a 
coupling function uj{4>) generalising the Brans-Dicke coupling constant [H]. This type of 
coupling is issued from particle physics theories whose Lagrangian at low energy could 
take the form of a scalar tensor theory. 

Geometrically, our present Universe seems well described by the isotropic and 
homogeneous cosmological models, i.e. the FLRW models. However the observations, 
as instance from Boomerang |31, show slight anisotropics in the cosmological microwave 
background, which could take origin at early times. Moreover, if the Universe had 
always been perfectly isotropic and homogeneous, it would be difficult to explain the 
large-scale structures we observe. Hence, it is interesting to consider an anisotropic 
Universe described by the Bianchi models. There are 9 ones but the most studied are 
the Bianchi type J, V, VIIq, Vllh and IX which are able to isotropise toward an 
FLRW model |n|. We will consider the Bianchi type / model which can tend toward 
a fiat FLRW one and is a good candidate from the infiation theory point of view. Of 
course the Bianchi models are not a definitive geometrical description of the Universe 
which should probably be inhomogeneous. But such models allow studying the necessary 
conditions for its isotropisation. 

Our goal will be to look for the necessary conditions depending on the potential and the 
Brans-Dicke coupling function for Universe isotropisation at late times. We will then 
derive the asymptotical dynamical behaviour of the metric functions and the condition 
for the presence of infiation. 

Technically, we will use the ADM Hamiltonian formalism [2llH| allowing to write the field 
equations as a first order system and then the dynamical systems theory as described in 



Isotropisation of Generalized Scalar-Tensor theory plus a massive scalar field. 



3 



[H] and suggested in [IS] to study them. We have not found any paper in the literature 
where these two methods are applied to equations system with 2 arbitrary functions. 
Due to this indeterminacy all the equilibrium points of the system can not be studied. 
However this problem can be overcomen for the subset of the phase space where lie the 
isotropic states of the Universe and which is of interest for us. 

This paper is organised as follow. In the second section we calculate the Hamiltonian 
field equations of the Generalized Scalar- Tensor theory plus a massive scalar field and 
rewrite them with new normalised variables. In the third section, we study the subset 
of the phase space corresponding to isotropy. We discuss about physical meaning of the 
mathematical results thus obtained in the fourth section. 



2. Field equations 

The Lagrangian of the Generalized Scalar- Tensor theory plus a massive scalar field is 
written: 

S = (IQn)-^ [\R- (3/2 + cj(0))0'^0>2 _ ^^d^x (1) 



with (j) the scalar field, u) the coupling function between the scalar field and the metric, 
U the potential. We will use the following form of the metric: 

ds^ = -{N^ - NiN')dn^ + 2NidVLj + RlgijJuj^ (2) 

the cj' being the 1-forms defining the Bianchi type / homogeneous space. The Qij are 
the metric functions, N and iVj are respectively the lapse and shifts functions. Using 
the methods described in [7] and fTT|, we find that the action can be rewritten in the 
following way: 

S = (167r)-i / {W^^ + n-^^ - NC'> - NiC')d^x (3) 

The n*-^ and II"^ are respectively the conjugate momentum of the metric functions and 
scalar field, the N and iVj play the role of Lagrange multipliers. The quantities Cq and 
Ci are respectively the super-Hamiltonian and supermomentum defined by: 

C = n« (5) 

where the hold for the quantities calculated on the 3-space and the "|" for the 
covariant derivative in the 3-space. When we vary the action with respect to the 
Lagrange multipliers we find the constraints C° = = 0. 

We rewrite the metric functions as Qij = e-2f^+2ftj g^j^^ ^j^g Misner 

parameterisation [8^: 

pl = 2nIV, - ^TTdlUl (6) 
6pij = diag{p+ + V3p^,p+ - VSp^, -2p+) (7) 
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Pij = diag{p+ + - V3(3., -2/5+) (8) 

Then, the action (jS} is written as: 

S = J p+dp+ + p.dp_ + p^d(f) - Hdn (9) 

with p^ = ttU^ and H = 27rn^. Finally, from the constraint C° = 0, we get the 
expression for the ADM Hamiltonian: 

H^=pI+ pi + 12^^ + 247r^R',e-^''U (10) 
3 -|- 2,uj 

from which we derive the Hamiltonian equations: 



dH 12(/)2 



dp^ (3 + 2iu)H 



(12) 



dH , , 

= -g^ = (13) 

d(P {3 + 2u)H (3 + 2cu)2if ^ H ^ ' 

A dot means a derivative with respect to fl. We will choose iV* = and we calculate N 
by writing that d^/dQ = -l/2n^A^ 0. Then, it comes: 

1 9TrR^p~^^ 

N = ^^""Y (16) 
H 

The relation between the VL and t times is then dt = —NdQ. We want to rewrite some 
of these equations under the form of an autonomous first order system with normalised 
variables p]. For this we will only consider the set of equations (fT2|l . (fT^ and (fT5|l . 
the equations (fTTni showing that the conjugate momentums of the variables describing 
the anisotropy are some constants. The constraint (fTn|) suggests the following set of 
normalised variables: 

X = H-^ (17) 

y = Ve-^^UH'^ (18) 

z = p^(P{3 + 2u)-^/^H-^ (19) 

The first one depends on if, the second one on H and (j) and the third one on H^ (p and 
P(j). Thus they are independent variables, y and z will be real if the functions U and 
3 + 2u; are positives. Under these conditions, it follows that the potential will favour 
infiation and the coupling function, when U = 0, will be such that the energy density 
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of the scalar field is positive. Rewriting the constraint (fTIHl with the new variables, we 
get: 



p'^x^ + K'y' + I2z 



2„.2 



(20) 



The positive constants p and are defined by p^ = p\+ p^_ and = 247r^i?Q. From 
this last equation we deduce that the variables x, y and z are bounded and belong to 
the following intervals: 

^ (21) 

(22) 



p ^,p ^ 



X G 
ye[-R,R] 



z e 



1/V12, 1/V12 



(23) 



The field equations (fT^ . (fTl|l and (fTKl) . become (see appendix [HI): 

X = 3R^y^x (24) 
y = y{Uz + 3i?2y2 _ 3^ (^25) 

i = y2/?2(3z-^£) (26) 

with (. = (pU^U^^^S + . They can not be expressed only with x, y and z because 

we do not wish to specify the form of uj and \J which are arbitrary functions of the scalar 
field. However, we do not need to know the exact form of £(x, y, z) since we are only 
interested by the asymptotical isotropisation of the Universe at late times. To reach our 
goal, it is sufficient to assume two types of asymptotical behaviours for i: ever it tends 
toward a constant or it diverges. This excludes any asymptotical chaotic behaviour for 
i and is in accordance with much of the functions u and U studied in the literature. 
We will have to check if these behaviours are compatible with the isotropisation of the 
Universe at late times. 



In the next section, we examine the equations ()24ll26p from the dynamical systems 
theory point of view. Firstly, we look for monotonic functions and secondly, we study 
the presence of equilibrium points. 



3. Dynamical studies of the fields equations 

Monotonic functions 



Lets examine the presence of monotonic functions. From the equation ^M . we deduce 
that a; is a monotonic function: when it is positive (negative), it increases (decreases). 
Since x has a constant sign, it follows from (fT8|l that it is the same for y. Notes also that 
in the plane x = with £ = etc, 2; is a monotonic and increasing function if z > £/6, 
decreasing otherwise. Thus there is no periodic or homoclinic orbit and then no chaotic 
behaviour. 



Isotropisation of Generalized Scalar-Tensor theory plus a massive scalar field. 



6 



If we look for the signs of the derivatives of the metric functions with respect to f2 
depending on the position of a point (x, y, z) in the phase space, we see that the sets 
of points such that they are constants are splat by planes defined by a; = cte since 
dgij/dVt = — 2e~^^"'"^^'J (1 — fiij). As instance for gu, it is defined by x = (p+ + \/?)p-)~^ 
and the sign of its derivative above or below this plane depends on the value of the con- 
stant p+ + v^p_. For (722, the plane is a; = (p+ — v^p_)~^ and for (733, x = p'^^ . Since 
X is a monotonic function with constant sign, we deduce that each metric functions can 
have one and only one extremum. From (|TK|l and the relation between Vt and t, we re- 
mark it will be the same in the proper time t. This is in agreement with the results of |Tn] 



Study of the isotropic equilibrium states 



Now we study the equilibrium points. They are all defined by {y,z) = (±(3 — 
£2^1/2^^2^2^-1/2 £^g^ and they will respect the constraint if x = 0. We have shown 
in [11]] that the Universe isotropises in the proper time t only when Q —00. This 
value of fl indicates that they will be sources or sinks but not saddle points. Thus an 
equilibrium states will represent an isotropic one for the Universe if in the same time 
f2 diverges negatively. Moreover, since x is a monotonic function of constant sign, we 
deduce from the relation ()16|) that the proper time is a monotonic and decreasing (in- 
creasing) function of fl when the Hamiltonian is positive (negative). Hence fl can be 
considered as a time variable and the equilibrium will take place at late times if H > 0. 
In what follows, we will assume that £ asymptotically tends toward a constant or di- 
verges. 

First, we assume that i is asymptotically a constant. We can show in that case by in- 
tegrating (I25II26I1 that when y = ±(3 - £^)^/^(3i?2)"^/^ Q diverges. It follows that these 
two equilibrium points are compatible with the isotropisation of the Universe. being 
a positive constant, they will be real if £^ < 3. We can not calculate their corresponding 
eigenvalues and thus knowing the signs of these last quantities because we do not now 
the expressions of the derivatives of i with respect to x, y and z. However, since we 
consider f2 as a time variable, they will be sinks if if > or sources if if < since they 
will respectively correspond to asymptotical late or early times. 

To get the behaviour of the metric functions when we approach an isotropic state, we 
need a differential equation for x when i — > cte. In this last case, the integration of the 
field equations p4ll26p gives: 



1 + GR'zo) - JP{1 + QR^Zo) + 18i?2zo(i?2y2 _ 1) 



(36i?2^o)"H27) 



By introducing this expression in the constraint equation and using (^Hl to express y 
as a function of x and its derivative, we get a differential equation for x. Since when 
the Universe isotropises, x and its derivative tend toward zero as VL diverges, and keep- 
ing only the second order terms in x and x, we find that x asymptotically behaves as 
XqCxp [(3 — £^)f2] when it tends to vanish. Taking into account the divergence of we 
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see that our approximation will be justified if P < 3, which is in accordance with our 
previous results. One could also recover this result by linearizing the equation ([2^ but 
we find this demonstration more rigorous. 

Now we examine the case for which £ diverges. It implies that the equilibrium points 
are unbounded. However, since y and z are bounded, we deduce that an equilibrium 
state can not be reached when i diverges. 

In the next section, we discuss about physical meaning of our results. 



4. Discussion 



In this work, we have examined the conditions under which the Universe described by 
a Generalized Scalar- Tensor theory with a massive scalar field isotropises as well as the 
asymptotical behaviour of the metric functions by help of the Hamiltonian formalism 
and dynamical systems theory. 

The set of points of the phase space corresponding to stable isotropic states for the Uni- 
verse is such that the time coordinate fl and the Hamiltonian diverges (Q ^ —oo and 
X — > 0). Then, the functions P± describing the anisotropy asymptotically tend toward 
a constant. We have shown that when i was asymptotically unbounded, an equilibrium 
state could not be reached. Thus, the isotropy of the Universe is not compatible with the 
divergence of the quantity 0[/<^?7~^(3 + 2u;)~^/^. However it arises when its value belongs 
to the range —a/3, \/S . In this case, the plane x = contains two equilibrium points 
corresponding to an isotropic state for the Universe. They are late times attractors in 
the t time if the Hamiltonian is a positive function. If H is interpreted as an energy, 
it means that we assume a positive energy for the Universe, which is reasonable. If it 
is not the case, the isotropisation arises at early times. Moreover, we have shown that 
as long as £^ < 3, the function x{fl) asymptotically tended toward Xoexp[{3 — 1"^)^]. 
Using Jini), we see that it corresponds to a power law of the proper times with the 
exponent if i does not tend toward a vanishing constant, or toward an exponential 
of t otherwise. These two types of functions represent the only possible attractors when 
the Universe isotropises. All this can be summarised in the following important result: 

A necessary condition for isotropisation of the Generalized Scalar-Tensor theory plus 
a massive scalar field (j), whatever the Brans-Dicke coupling function oo and the potential 
U considered, will be that (pU^U^^ {3 + 2uj)^^/'^ tends toward a constant i with P' < 3. It 
arises at late times if the Hamiltonian is positive, at early times otherwise. If £ ^ the 
metric functions tend toward t^ ^ . The Universe is expanding and will be inflationary if 
(? <1. If £ = 0, the Universe tends toward a De Sitter model. 



Note, that the asymptotical behaviour of the metric functions when isotropisation arises 
does not depend on initial conditions whereas the epoch of isotropisation, i.e. late or 
early times, depends on the initial sign of the Hamiltonian. One element is missing 
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in this result: the value of the scalar field when the isotropisation arises, i.e. when 
Q —oo. Expressing as a function of z and (see appendix EJ, and taking z as its 
value at the equilibrium, £/6,we get a differential equation for 0. It does not describe 
the scalar field behaviour during the whole Universe evolution, but asymptotically when 
fl —oo and the system approach equilibrium. This equation of the first order can 
be solved analytically or numerically. This additional result allows to calculate i when 
isotropisation occurs and completes the main one above. It is written: 

The value of the scalar field when the Universe reaches an isotropic equilibrium state is 
the value of the function defined by (p = 2(f)^U^{?> + 2uj)^^U^^ when VL —oo. 

Lets examine the relations between these results and others quoted in the literature. 

Firstly they are in accordance with the "No Hair Theorem" which states that 
General Relativity with a cosmological constant isotropises toward a De Sitter model 
since in this case, £ = 0. It will be true for any form of potential and Brans-Dicke 
coupling function such that ^ asymptotically vanishes when Vt —>■ — oo which does not 
necessary implies that the potential tends toward a constant. As instance, it arises if 
the Brans-Dicke coupling function diverges faster than (pU^U'^ . This generalises the 
"No Hair Theorem" in the special case of the Bianchi type / model. 

Secondly, in ^2], it has been shown that all the Bianchi models with an exponential 
potential V = e'"^ (except the contracting Bianchi type IX model), isotropise at late 
times when k"^ < 2. If = 0, these models tend toward a De Sitter model and toward 
t^'^ ^ otherwise. If /c^ > 2, the Bianchi type J, V VII and IX models might isotropise 
at late times. In the present paper, the form of the coupling constant corresponding to 
the theory studied in [12] is ^/3~+~2cJ0^^ = \/2. What can we deduce from our results? 
If we introduce these forms of u and U in the expression of we see that the necessary 
condition for the isotropisation of the Bianchi type / model will be k"^ < 6. Then, the 
Universe is of De Sitter type if = 0. In the other cases, the metric functions behave 
as t^^ ^ The inflation arises when fc^ < 2. All these results are in accordance with 
these of [12] and [12]. However, some differences exist which are not in contradiction 
with the previous quoted papers: we have shown that Universe might isotropise and is 
inflationary when k"^ <2 but isotropisation is impossible if A;^ > 6. Between these two 
values, the necessary condition for isotropy is respected but no inflation can occur. 

Last, these results are agreed with these found in [H]. In this paper where the 
Hyperextended Scalar Tensor theory with a potential is studied for the Bianchi type I 
model, it is shown that the Universe isotropises when / Ge^^dt tends toward a constant, 
G being the gravitational coupling function. If in this last expression we choose G = 1 
and e^^ t~^^ ^, we find that isotropisation arises if < 3, in agreement with the 
above results. 

Lets say few words about the power law potential, U = (p^. We can show from the 
asymptotical equation for defined above, that when Q —oo, +oo if A; < (if 
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A; > 0, the scalar field is not real). Then £ ^ and isotropisation systematically leads 
to an asymptotical De Sitter model. 

The result of this paper is not only a necessary condition for the isotropisation of 
the Universe. We have also derived the asymptotical behaviour of the metric functions 
and thus some conditions for a late time inflationary behaviour. It is a strong theoretical 
constraint on the forms of u and U so that the Universe be physically realistic at late 
times if we consider that it can be described by a Generalized Scalar- Tensor theory with 
a massive scalar held in the Bianchi type / model. We have checked the compatibility 
of our results with these of the important No Hair Theorem and these of Kitada et al 
and Cooley et al that are here unified in a single condition. To our knowledge, there 
is no paper mixing the Hamiltonian technique and the dynamical systems theory with 
so many arbitrary functions. It seems to be a fruitful method in the case studied here 
mainly because it allows to calculate the equilibrium points as function of the potential 
and Brans-Dicke function. Then from mathematical constraints on the equilibrium 
points, we deduce constraints for these undetermined quantities. In future papers, we 
will see that we get the same type of results when we introduce a perfect fluid and we 
will extend this method to more general theory such has Hyperextended Scalar Tensor 
ones or other Bianchi models. 

5. Appendix 

Equation for x 

From (|T5| and (fTTjl . we deduce: 

X = SR^y^x (28) 

Equation for y 

By deriving (fT8|l . we find: 

<i) = 2L^'x~^y{y-3RV) (29) 
with L2 = {e~^^U)(f,. By using (fT2l) and (fT9ll to express p^, we find: 

= 1 (30) 

with L3 = (3 + 2ujy/^(t)-^ Consequently, ^ and JHOI) gives: 

y = 6L2L^^x^y-^z + 3i?V (31) 
Equation for z 

by using (fT9|) to express and by deriving this expression, we get: 

p^ = L^x-H - SR^L^x-^y^z - U^'^x'^z^ + l2LiLf<p~'^x'^ z^ (32) 
with L4 = uj^. From (fT^ and by using the fact that U~^U^ = L2x'^y''^ + 1/2^32;"^, it 



comes: 



p^ = -I2(j)-^x~^z^ + l2ULl'^<p-'^x~^z^ - -7^L2X - —Lsx'^y^z-^ (33) 
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From the equations (p^ and (jSBI), we derive: 

i = ^RV^ - ^L,L^'x' - ^y'z-' (34) 

Before getting the equation IHM to we have to evaluate the term L-zL'^^x^y''^ . 

After few calculations, we find -{2z)'^ + 0f/-if/<^(3 + 2uj)-^''^. 
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